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ABSTRACT 


Au  analysis  based  primarily  on  computations  is 
conducted  to  determine  errors  in  a bombing  system 
similar  to  that  currently  under  test  and  development 
by  the  Marine  Corps. 

There  are  many  parameters  in  the  system 
determining  point  of  impact.  Included  are  point  of 
release  and  velocity  at  that  point,  air  density, the 
mass  and  drag  coefficient  of  the  bomb,  its  cross 
sectional  area,  and  the  wind  structure.  Most  values 
are  not  known  with  the  desired  accuracy  but  must  be 
estimated,  say,  from  radar  measurements,  for  position. 

This  study  estimates  as  much  as  possible  the  error 
in  predicting  the  point  of  impact  as  a consequence  of 
such  factors  as  random  uncorrelated  noise  in  the  radar 
measurements,  errors  in  estimating  density,  errors  in 
the  wind  structure,  and  errors  due  to  initial 
conditions. 


4 


TABLE  OF  CONTENTS 


1.  INTEOOUCTION 7 

II.  EfiBOfiS  DUE  TO  CHANGES  IN  DENSITY,  NINE  STBOCTURE,  AND 
VALUE  OF  DBAG  COEFFICIENT 10 

A.  EQUATIONS  OF  MOTION  FOB  A PROJECTILE 11 

B.  BASIC  ASSUMPTIONS  AND  CONDITIONS 12 

1.  Altitude  and  Airspeed 12 

2.  Drag  Coefficient  and  Constant  Values 12 

3.  Air  Density 12 

4.  Hind  Structure 13 

C.  COMPARISON  OF  SEVERAL  HIND  CONDITIONS 13 

D.  COMPARISON  OF  VALUES  OF  THE  CONSTANT  C 20 

III.  EFFECTS  OF  INITIAL  CONDITIONS  ON  POINT  OF  FALL...  22 

A.  VARIATIONAL  EQUATIONS 23 

1.  Effects  of  Changes  in  v 25 

2.  Effects  of  Other  Changes 25 

B.  THE  ADJOINT  SYSTEM 28 

1.  Derivation  of  Influence  Coefficients 26 

2.  Analysis  of  Results 31 

C.  COMMENTS 34 

IV.  ANALYSIS  OF  FACTORS  INFLUENCING  AIRCRAFT  POSITION  AND 

VELOCITY 35 

A.  RANDOM  UNCORRELATED  NOISE 35 

1.  A Nunerical  Approach 36 

2.  An  Analytic  Approach 36 

3.  Comparison  of  Results  Obtained 41 

B.  PERIODIC  ERRORS 41 

1.  The  Analysis 41 

2.  The  Results 44 

C.  EFFECTS  OF  TURNS 46 

V.  CONCLUSIONS  AND  COMMENTS 48 


5 


A.  CONCLUSIONS 46 

B.  COMilENXS 50 

LIST  OF  BEFEBENCBS 52 

INITIAL  DISTRIBUTION  LIST 53 


6 


I,  INTRODUCTION 


This  thesis  developes  and  uses  oany  of  the  mathematical 
techniques  and  relations  that  are  needed  to  analyze  the 
important  errors  in  a bombing  system.  The  AN/TPQ’-27, 
described  in  Ref.  1 and  Ref.  3,  is  the  model  for  the 
analysis. 

Precision  bombing  is  a crucial  means  of  supporting 
Marine  Corps  operations.  It  is  a matter  of  life  and  death 
that  it  be  as  accurate  and  effective  as  possible.  To  ensure 
this  the  Marine  Corps  carries  out  extensive  analysis  and 
testing  of  its  systems. 

The  complete  analysis  consists  of  two  basic  parts. 
First  there  is  a mathematical  modeling  of  the  system;  that 
is,  equations  must  be  formulated  to  describe  its  behavior 
and  interreactions  of  its  components.  These  must  be  based 
on  the  natural  laws  that  govern  the  motions  and  on  the 
accumulated  experiences  (some  lore  perhaps)  with  similar 
systems.  The  equations  must  be  simple  enough  to  allow  the 
necessary  solutions  and  computations,  but  they  must  also 
maintain  the  needed  accuracy.  This  thesis  takes  up  this 
part.  Second,  there  must  be  extensive  testing  of  the 
components  and  of  the  integrated  system  under  conditions 
which  simulate  the  field  as  closely  as  possible,  but  with 
extra  measurements  made  to  isolate  and  determine  the  sources 
of  significant  causes  of  error. 

In  this  study  the  author  nas  set  up  the  basic 
mathematics  for  error  analysis  and  carried  out  computations 
to  show  the  typical  effects  of  errors  caused  or  observed  at 
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various  stages 


It  is  convenient  to  consider  three  classes  of  errors, 
depending  on  the  time  they  are  observed.  Different  analysis 
is  reguired  for  each  type. 

First  are  the  errors  which  are  due  to  inaccuracies  in 
the  equations  of  motion  during  the  period  of  fall.  These 
may  come  from  sources  such  as'  errors  in  estimating  the  drag 
coefficient,  and  variations  in  the  air  density  and  the  wind 
structure  from  the  assumed  values. 

Second,  there  may  be  significant  errors  in  the  position 
and  velocity  at  the  time  of  release.  These  may  come  from 
sources  such  as  wind  gusts  near  the  time  of  release,  faulty 
plane  response,  errors  in  the  filtered  estimates  of  position 
and  velocity,  and  release  or  separation  problems. 

Third,  there  may  be  errors  in  the  radar  measurements  and 
transients  introduced  in  various  ways.  The  errors  may  be 
inherent  in  the  radar,  they  may  be  due  to  the  changing 
aspect  of  the  plane,  to  maneuvers,  and/or  to  the  filtering 
equations  whose  function  is  to  reduce  the  effects  of  random 
errors.  We  will  know  only  statistical  properties  of  some  of 
these,  and  for  these  we  can  only  determine  associated 
statistical  properties  of  the  associated  miss  distances. 

A number  of  estimates  were  made  in  two  or  more  ways,  for 
comparison  of  the  methods  and  to  act  as  a check.  For 
example  the  effects  of  variations  in  initial  velocity  were 
estimated  in  three  ways:  (1)  by  recalculating  the 
trajectory,  (2)  by  using  variational  equations,  and  (3)  by 
using  the  adjoint  or  costate  equations.  For  small  changes 
the  latter  two  are  more  accurate;  for  large  changes,  the 
former  is  better,  but  all  are  adequate  over  a considerable 
range.  Estimates  of  errors  in  the  filtered  values  of 
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position  and  velocity  were  also  obtained  in  two  basic  ways: 
(1)  by  an  analysis  of  the  solutions  to  the  difference 
equations  that  are  used  for  filtering,  and  (2)  by  a 
numerical  simulation  of  the  solution,  on  the  computer. 

He  need  to  make  some  remarks  about  notation  and 
coordinate  sets  at  this  time.  He  will  use  an  z-y-z 
Cartesian  coordinate  set  fixed  with  respect  to  the  earth. 
The  program  has  a subroutine  which  corrects  automatically 
for  the  eeirth's  rotation  so  we  can  consider  the  set  to  be 
Newtonian.  The  z direction  is  the  vertical.  The  plane  is 
assumed  to  be  flying  horizontally,  and  its  velocity  (with 
respect  to  the  ground)  defines  the  y-azis. 

Derivatives  with^respect  to  time  are  indicated  by  a dot 
over  a variable,  X * , etc.  In  the  tezt,  however, 
this  is  denoted  by  zdot,  and  a similar  notation  is  used  for 
other  variables. 

The  ratio  is  used  in  the  equations  as  the  ratio  of 
the  air  density  at  altitude  z to  the  air  density  at  ground 
level  (He  use  an  exponential  value  for  this.).  In  the  text 
these  terms  are  written  rho  and  rho(O)  respectively.  In  the 
same  way  variations  such  as  are  written  as  DELx  in  the 
text. 

It  was  not  feasible  to  write  subscripts,  so  that 
subscripted  variables  such  as  the  drag  coefficient,  the  z 
component  of  the  wind  velocity,  the  y component  of  the  wind 
velocity,  and  others,  were  written  Cd,  Vwx,  Vwy,  etc.,  in 
the  tezt;  it  is  felt  that  this  will  not  be  confusing. 

The  notation  was  chosen  to  be  in  conformity  with  that  in 
Bef.  1 and  Ref.  3 as  much  as  possible.  Also,  the  terms 
impact  and  fall  are  both  used  interchangeably  to  denote 
conditions  at  the  instant  the  bomb  strikes  the  ground. 
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II.  ERRORS  DUE  TO  CHANGES  IN  DEN SI lY.  HIND  STRUCTURE.  AND 
VALUE  OF  DMG  COEFFICIENT 


Hhen  the  ordnance  is  released  it  follows  a path 
deterained  by  natural  laws.  He  approximate  these  as  best  we 
can  to  predict  the  point  of  impact.  Houeverf  the  equations 
of  motion  contain  certain  parameters  and  functions  which  may 
not  be  known,  and  the  planned  trajectory  has  a point  of 
impact  different  from  that  calculated,  and  a miss  results. 

The  effects  of  two  of  these  parameters  are  investigated 
here.  The  first  of  these  is  the  wind.  Let  us  assume  that 
we  have  a good  estimate  of  the  wind  at  the  point  of  release 
(called  the  wind  at  altitude) , because  we  know  the  velocity 
of  the  plane  (from  the  radar) , its  airspeed,  and  its 
heading.  However,  the  wind  varies  with  altitude,  and  the 
structure  of  this  function  may  be  unknown  under  combat 
conditions.  Hhen  testing,  we  can  approximate  this  function 
by  making  extra  meteorological  measurements.  Then  by  using 
various  approximations  to  the  wind,  based  on  the  wind  at 
altitude,  we  can  calculate  the  associated  miss  distance.  An 
approximation  that  is  commonly  used  is  that  the  wind  is 
constant  or  uniform,  that  constant  being  some  fraction  of 
the  wind  at  altitude.  The  lore  indicates  that  eight-tenths 
is  a good  fraction;  higher  values,  approaching  one,  have 
been  suggested  [Ref.  4].  Under  the  conditions  of  our 
example  eight-tenths  seems  like  a good  choice.  However,  the 
component  of  the  wind  normal  to  the  wind  at  altitude  often 
causes  a large  error  that  there  is  no  way  to  correct  knowing 
just  the  wind  at  altitude.  Only  a further  study  of  wind 
structures  will  show  when  bombing  is  ineffectual  without 
more  detailed  knowledge  of  the  wind. 
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A second  factor  whose  effects  we  can  estimate  readily  is 
the  drag  parameter.  There  is  a single  drag  parameter,  which 
involves  the  air  density  at  some  reference  level,  the  cross 
sectional  area  of  the  bomb  and  its  mass,  and  a drag 
coefficient,  Cd,  aepending  primarily  on  its  shape.  Ne  can 
estimate  the  effects  of  errors  in  this  parameter  directly. 

There  are  many  other  sources  of  error  for  which  no 
simple  analysis  exists;  some  of  these  are  discussed  later  in 
this  chapter. 

A.  EQUATIONS  OF  ^lOTION  FOR  A PROJECTILE 


Reference  1 describes  in  detail  the  equations  involved 
in  the  ballistic  integration.  These  are  the  classical 
equations  of  motion  for  a non-lifting  projectile: 

X-  -/l/~^c^l/Cx-yJ  ('-> 

C^V  (f-V^)  « 

-j  - Aj^x  ^ !/  4 Oc-) 


where  A is  the  cross  sectional  area  of  the  projectile,  m is 
the  mass  of  the  projectile,  g is  the  acceleration  due  to 
gravity;  rho,  Cd,  xdot,  ydot,  zdot,  Vwx,  and  Vwy  have  been 
previously  defined;  and  V,  the  airspeed,  is: 


Equations  (la),  (1b),  and  (1c)  were  integrated 

numerically  to  obtain  final  x,  y,  and  z coordinates  and 

* 
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hence,  to  estimate  the  bomb  impact  point 


B.  B&SIC  ASSUMPTIONS  AND  CONDITIONS 


It  vas  decided  to  use  values  for  various  parameters  that 
were  selected  arbitrarily  in  a reasonable  range,  rather  than 
treat  a specific  missle,  whose  details  would  be  classified. 

Altitude  and  Airspeed 

In  all  cases  the  aircraft  was  assumed  to  be  flying 
straight  and  level  at  an  altitude  of  20,000  feet  with  an 
airspeed  of  350  feet  per  second.  In  those  cases  where  a 
wind  structure  was  considered  the  initial  velocity  (with 
respect  to  the  earth)  was  altered  to  take  into  account  the 
wind  velocity. 


2.  Drag  Coefficient  and  Constant  Values 


The  drag  coefficient,  Cd;  the  cross  sectional  area 
of  the  projectile.  A;  the  mass  of  the  projectile,  m;  and  the 
air  density  at  ground  level,  rho(O) ; were  assumed  for  the 
majority  of  the  work  to  take  on  values  such  that  the 
important  ratio  C=(A/m)  (Cd)  (rho(O))  had  the  value  0.000025. 
This  figure  was  considered  to  be  a good  approximation  to 
some  typipal  ordnance  types  without  requiring  the  thesis  to 
be  classified. 

3.  iir 

It  was  assumed  the  decrease  in  air  density  with 
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altitude  was  given  by  the  factor  exp(-z  ln2/18000) . This 
function  was  compared  with  the  ratio  given  in  Ref.  1,  and  it 
was  felt  that  the  exponential  expression  was  satisfactory 
(a:>d  simplified  programming)  . 

4.  Hind  Structure 


Reference  4 details  wind  structure  data  obtained  by 
Camp  Pendleton  Harine  Corps  Base  in  1975.  The  data  was 
collected  twice  daily  on  25  March,  28  March,  and  12  April. 
Hind  velocity  in  knots  and  wind  angle  were  determined  at 
2,000  foot  intervals  from  ground  level  to  36,000  feet.  All 
reference  within  this  work  to  a "variable  wind  structure" 
applies  to  the  data  in  Ref.  4 collected  25  March,  1975,  at 
11:00  unless  otherwise  noted. 


C.  COMPARISON  OF  SEVERAL  HIND  CONDITIONS 


Several  trajectories  were  computed  based  on  two  actual 
wind  structures.  In  each  case  directions  were  chosen  for 
the  plane  corresponding  to  a headwind,  a crosswind,  and  a 
tailwind.  In  each  case  trajectories  were  computed  with  a 
uniform  wind  equal  to  the  wind  at  altitude,  and  one  equal  to 
eight-tenths  wind  at  altitude.  In  two  cases  various  other 
uniform  winds  were  also  assumed. 


and 


X 


In  all  calculated  trajectories  the  following  equations 
terms  were  used: 


^ Urn 
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w 

(3c) 

l/=  ^ 

CsJ) 

o.  ooooz^ 

O*-) 

• po  e,x.p(-^  (^z//rff0o)^ 

Of) 

^■■CfV(..lc,) 

O3) 

v.-CfV(.-y^^) 

Oi) 

~ ~ Cj>  V w 

00 
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All  trajectories  were  calculated  using  a modified  Euler 
Method  £Ref.2]  with  a time  step,  DI=1.0. 

Since  the  wind  velocity  components  in  Ref.  4 are 
provided  at  2,000  feet  intervals  only,  a linear 
interpolation  scheme  was  devised  to  compute  wind  velocity 
components  as  the  bomb  fell. 

Three  trajectories  using  the  wind  structure  at  11:00,  25 
March,  1975,  were  then  computed  for  three  conditions:  a 
plane  flying  into  a headwind,  a plane  flying  wixh  a 
crosswind,  and  a plane  flying  with  a tailwind.  Next,  wind 
at  altitude  was  considered  as  the  constant  wind  velocity 
throughout  as  opposed  to  using  a variable  wind  structure, 
and  the  three  trajectories  computed  again.  Finally,  the 
three  trajectories  were  computed  again  assuming  the  wind  to 
be  uniformly  eight-tenths  that  at  altitude. 

The  trajectories  based  on  the  variable  wind  structure 
ware  used  as  trajectories  of  reference  and  the  errors  due  to 
the  approximations  calculated.  Corresponding  trajectories 
using  the  stratified  wind  data  of  25  March,  1975,  14:00, 
contained  in  Ref.  4 were  then  calculated  under  the 
conditions  described  above,  for  a comparison. 

If  the  wind  at  altitude  is  icnown,  and  is  aj.1  -that  is 
known  about  the  wind,  the  only  uniform  wind  we  may 
reasonably  choose  is  one  which'  is  some  fraction  of  the  wind 
at  altitude.  It  was  felt  that  the  change  in  impact  point 
due  to  changing  this  fraction  would  be  approximately  linear, 
lo  check  this  two  trajectories  were  calculated  using 
nine-tenths  and  zero  respectively,  as  fractions.  The 
resulting  trajectory  using  nine-tenths  led  to  a point  of 
impact  which  was  approximately  the  midpcint  between  the 
points  due  tc  fractions  of  eight-tenths  and  one,  indicating 
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the  change  was  appcoximately  linear.  A typical  plot  (Figure 
1)  showed  that  in  the  example  considered,  the  value  of 
eight-tenths  was  near  the  optimum.  It  is  felt  no  other 
value  would  result  in  a significantly  smaller  error  in  any 
of  the  cases  considered. 

While  Sef.  4 points  out  that  fractions  of  wind  at 
altitude  closer  to  unity  are  better  for  certain  bomb  loads, 
the  results  indicate  that  for  these  wind  structures, 
eight-tenths  is  near  the  best  choice  for  use  in  the  absence 
of  stratified  wind  data.  The  resulting  error  was  near  the 
minimum  in  all  cases  for  this  choice. 

For  these  wind  structures  there  was  always  a sizeable 
error  when  a constant  wind  was  assumed,  due  to  changes  in 
the  direction  of  the  wind  with  altitude.  Hore  analysis  with 
different  wind  observations  is  needed  to  establish  the 
validity  of  such  a universal  factor  and  the  general  errors 
due  to  the  varying  direction  of  the  wind. 

Some  additional  results  of  the  analysis  are  contained  in 
TABLES  1 and  II. 

In  the  next  chapter  another  method  is  given  for 
estimating  the  effects  of  the  changes  discussed  above,  using 
methods  like  those  developed  by  G.  A.  Bliss  for  corrections 
in  artille'ry. 
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Figure  1 


A typical  Plot  of  the  Errors  From  the  Variable 
Trajectory  Impact  Point  (in  Ft.)  When  Various 
Fractions  of.  Wind  at  Altitude  Are  Assumed 
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oTRATIPIEU 

WIND  data 

FROM  25  March. 

75.  11: 

00 

type 

wind 

wind  vel 
used 

X 

y 

X error 

y error 

tail 

variable 

67.754 

15519.350 

N/A 

N/A 

tail 

wind  at 
alt 

-.028 

15629.180 

-67.8 

1 10 

tail 

.8  wind 
at  alt 

-.023 

15499.350 

-67.8 

-20 

cross 

variable 

-409.653 

10994.860 

n/a 

N/a 

cross 

wind  at 
alt 

-518.572 

10922.150  - 

108.9 

-72.1 

cross 

.8  wind 
at  alt 

-412.139 

10929.470 

-2.5 

-65.4 

head 

variable 

-65.814 

7488.727 

N/A 

n/a 

head 

wind  at 
alt 

0.0 

7376.344 

65.8 

-112.4 

head 

.8  wind 
at  alt 

0.0 

7507.988 

65.8 

19.3 

Table  I 

X and  Y Errors  (in  Ft.)  for  Wind  at  Altitude 
and  .8  Wind  at  Altitude  Trajectories  Versus 
Variable  Trajectory  Under  Various  Wind  Conditions 
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STRATIFIED 

WIND  Data 

FROM  25  MARCH. 

75, 

00 

type 

wind 

wind  vel 
used 

X 

y 

X error 

y error 

tail 

variable 

57.813 

15865.250 

N/a 

N/A 

tail 

wind  at 
alt 

-.028 

16010.070 

-57.8 

144.8 

tail 

.8  wind 
at  alt 

-.023 

15868.360 

-57.8 

3.11 

cross 

variable 

-422.153 

10874.430 

N/A 

N/a 

cross 

wind  at 
alt 

-566.438 

10806.460  - 

144.3 

-68.0 

cross 

.8  wind 
at  alt 

-449.604 

10815.100 

-27.5 

-59.0 

head 

variable 

-56.073 

7142.984 

N/a 

n/a 

head 

wind  at 
alt 

0.0 

6995.461 

56.0 

-148.0 

head 

.8  wind 
at  alt 

0.0 

7139.324 

56.0 

-3.7 

Table  II 

X and  Y Errors  (in  Ft.)  for  Wind  at  Altitude 
and  .8  Wind  at  Altitude  Trajectories  Versus 
Variable  Trajectory  Under  Various  Wind  Conditions 
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0.  COHPABISON  OF  VALUES  OF  THE  CONSTANT  C 


The  value  C=  (A/m)  (Cd)  (rho  (0) ) =0. 000025  was,  as  stated 
earlier,  used  for  the  majority  of  the  computation.  In  this 
section  this  parameter  was  varied  to  determine  the  effects. 

Trajectories  were  computed  with  vario-us  values  of  C, 
namely  0-,  .0000125,  .000025,  and  .000050,  to  determine  the 
associated  change  in  point  of  impact. 

As  indicated  in  TABLE  III,  the  function  is  nearly  linear 
over  a wide  range,  and  hence  changes  in  point  of  impact 
could  be  estimated  by  linear  approximations. 
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c 

Y 

CHANGE  IN 
IMPACT  POINT 

.000025 

11502.730 

n/a 

.000050 

10707.760 

794.97 

.0000125 

1 192S.490 

425.7.6 

.0 

12373.250 

870.5 

Table  III 

Change  in  Point  of  Impact  (in  Ft.)  as  Result 
of  Varying  the  Constant  C (Trajectory  with 
C = 0.000025  is  Nominal  Trajectory) 
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III.  EFgECTS  OF  INITIAL  CONDITIONS  ON  POINT  OF  FALL 


Fron  the  point  of  view  of  control  the  time  of  release  of 
the  bomb  is  the  crucial  instant.  Conditions  at  this  time 
determine  success  or  failure,  and  nothing  further  can  be 
done  to  improve  bombing  accuracy. 

In  this  chapter  a study  is  made  of  the  changes,  or 
errors  in  the  point  of  fall,  as  a consequence  of  changes  or 
errors  in  initial  conditions.  By  initial  conditions  are 
meant  the  conditions  at  the  instant  of  release  from  the 
plane. 

The  effects  of  changes  in  the  initial  values  of  x and  y 
on  the  point  of  impact  are  obvious.  The  effects  of  changes 
in  u,  V,  w,  and  z are  not  so  obvious;  these  effects  are 
estimated  in  three  ways  and  the  results  compared.  The  first 
way  is  to  take  a typical  trajectory  and  then  recompute  the 
trajectory  several  times,  each  time  changing  one  initial 
condition  by  a "small  amount".  The  second  way  is  to  use  the 
variational  equations.  These  ware  integrated,  using  the 
appropriate  initial  conditions  in  order  to  derive  tha 
differential  estimates  of  the  changes  in  the  point  of 
impact.  The  third  way  is  to  define  the  adjoint  or  costate 
equations.  These  are  integrated  backwards,  with  appropriate 
terminal  conditions,  to  define  the  so-called  influence 
coefficients,  which  may  be  interpreted  as  the  partial 
derivatives  of  the  final  conditions  with  respect  to  the 
initial  conditions. 

Comparisons  of  the  results  were  carried  out  with  three 
reference  or  nominal  trajectories.  In  one  trajectory  there 
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was  assuaed  to  be  no  wind;  in  the  second  the  wind  was 
assumed  to  be  uniformly  eight-tenths  that  at  altitude;  in 
the  third  the  wind  function  was  assumed  to  be  that  observed 
at  Camp  Pendleton  on  25  March,  1975,  at  11:00. 

In  all  cases  all  of  the  methods  gave  very  similar 
results. 

A.  VARIATIONAL  EQUATIONS 

The  variational  equations  are  derived  as  follows. 

Consider  the  equations  of  motion: 

C.P  1/ 1/  } 

‘ 1>/X.  ' 

-C/>V(y-l/  ) (">) 

UJ  - - Cj)  V UJ  - ^ (ic) 

where  u,  v,  w,  C,  rho,  V,Vwx,  and  Vwy  have  been  previously 
defined.  Then  let  us  consider  the  first-order  effects  if  we 
change  the  variables  by  small  amounts  in  some  way,  but 
require  that  the  new  variables  satisfy  the  equations  of 


motion.  The  first-order  terms  in  the  change, 
must  then  satisfy  the  equations: 

OELz,  etc.. 

J'X  ’^fU. 

a^) 

M) 

Cecj 
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These  are  called'  the  variational  equations.  It  nay  be 
noted  that  they  are,iinear  but  with  variable  coefficients 
whose  values  oust  be  obtained  from  the  original  trajectory. 
If  the  variations  In  the  wind  are  zero,  the  equations  are 
also  homogeneous. 


1 . Effects  of  Changes  in  v 


As  a first  example  let  us  consider  the  effects  of 
changes  in  the  initial  speed,  v (0)  . Let  us  taice  as  a 
standard  or  nominal  path  one  with  no  wind.  The  values  of  z, 
u,  V,  w when  they  occur  in  the  variational  equations  are 
from  this  path. 

Let  us  consider  DELv(0)  = 1,  and  all  other  variations 
to  be  zero  initially.  We  will  integrate  to  the  terminal 
time  tf,  the  time  of  fall  on  the  original  path.  We  find 
that  DELy  (tf ) -31 . 30.  However,  DELz (tf ) = 1 . 34,  and  hence  the 
differential  change  in  terminal  time  (denoted  Dtf)  is 
Otf =”DEL2  (tf ) /w (tf ) =0 . 001 . The  differential  change  in  y 
(denoted  Dyf)  then  requires  a secondary  term  due  to  the 
change  in  terminal  time;  Dyf =DELy (tf) ♦v (t f ) (Dtf ) = 31 . 63.  We 
see  that  the  change  in  terminal  time  has  a secondary  effect. 

The  variation  in  x in  this  case  is  zero. 


2*  Effects  of  Other  Changes 


In  a similar  manner  the  effects  of  changes  in  other 
initial  conditions  were  studied.  The  variations  DElu(O), 
OELz(O),  and  DELw(O)  were  each  separately  set  equal  to  one 
initially,  while  all  of  the  other  variations  were  set  to 
zero  initially;  we  integrated  to  terminal  time,  tf,  and 
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computed  Dxf  and  Oyf  in  each  case.  We  then  computed  the 
error  in  impact  point  due  to  the  differential  changes  in  x 
and  y. 


Identical  calculations  were  performed  considering 
the  three  separate  nominal  trajectories,  and  as  indicated  in 
TABLE  IV,  differences  in  figures  among  the  three 
trajectories  are  quite  minimal. 
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NOMINAL 
■i’RA  J . 

iiRROK  V/HEN 
DELV(0)= 

1.0 

error  when 
v(0) 

INCREASED 

10.0 

ERROR  WHEN 
DELZ(0)= 

1.0 

ERROR  WHEN 

z(o) 

INCREASED 

10.0 

no  wind 

31.633 

315.938 

.283 

2.85 

.8  wind 
at  alt 

31.767 

317.0 

.246 

2.495  . 

variable 

wind 

31.671 

316.31 

.245 

3.845 

ERROR  WHEN 
DELW(0)= 

1.0 

ERROR  WHEN 
w(o) 

INCREASED 

10.0 

ERROR  WHEN 
DELU ( 0 ) = 

1.0 

ERROR  WHEN 

u(o) 

INCREASED 

10.0 

no  wind 

9.535 

96.305 

32.867 

328.646 

.8  wind 
at  alt 

8.372 

84.532 

32.799 

328.0 

variable 

8.306 

83.60 

32.728 

327.4 

wind 


Table  IV 

Comparison  of  Errors  in  Impact  Point 
(in  Ft.)  When  Errors  in  Various 
Initial  Conditions  are  Assumed 
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B.  TH£  ADJOINT  SYSTEN 


The  adjoint  system  Eor  the  equations  of  motion  vas 
derived,  and  the  influence  coefficients  analyzed  as  a means 
of  cross  checking  the  results  of  the  variational  equation 
analysis.  The  nominal  trajectory  for  this  comparison  had  no 
wind. 


1.  Derivation  of  Influence  Coefficients 


Let  us  introduce  six  new  unspecified  variables 
P 1 ,P2,P3 , P5,and  B6  to  be  defined  later.  The  hamiltonian 
for  the  system  is  defined: 

hi  - UL  + -f- - f>  C\/ 

where  u,  v,  w,  rho,  C,  Vwx,  and  Vwy  have  been  previously 
defined.  Now,  let  us  always  choose  the  P's  to  be  solutions 
to  the  differential  equations: 


where  xl=x,  x2=y,  x3=z,  x4=u,  x5=v,  and  x6=w. 


Then 

o.  zdiL 

' ' 3ic 


o 
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This  does  not  define  a solution  yet.  If  we  define 
values  for  all  six  components  at  some  particular  time,  a 
particular  solution  will  be  determined  completely,  though  we 
may  have  to  integrate  numerically  to  obtain  it. 


The  adjoint  variables  and  the  variations  satisfy  an 
important  identity 


which  is  easily  verified,  whenever  the  DELxi  are  solutions 
to  the  homogeneous  variational  equations. 


This  relation  is  the  basis  for  obtaining  what  are 
sometimes  called  influence  coefficients  or  sensitivity 
factors. 


To  illustrate  how  these  are  obtained,  let  us  derive 
the  sensitivity  factors  for  the  final  value  of  y,  or  x2. 
First,  let  us  observe  that  if  we  integrate  equation  (6)  from 
zero  to  tf  we  get 


Now,  let  us  choose  a solution  to  the  adjoint  such  that  at 
time  tf  the  coefficient  of  DELy=D£Lx2  is  P2(tf)=1,  but  all 
other  Pi  (tf ) are  zero.  To  generate  this  solution  we  must 
integrate  backwards  to  t»0.  For  this  particular  solution  to 
the  adjoint  system  we  get  then  from  equation  (7) 
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That  is,  for  this  particular  solution  the  components  Pi(0) 
yield  the  influence  of  DELxi(O)  on  DELy(tf),  and  hence  the 
name. 


Several  comments  may  be  made.  There  is  one  such 
solution  vector  Pi  for  each  variable  z,  y,  z,  u,  v,  and  w. 
These  were  calculated  for  their  corresponding  variables. 
The  values  obtained  in  this  way  for  the  final  variations 
should  be  the  same  as  those  obtained  from  the  variational 
equations  except  for  errors  in  numerical  integration  due  to 
roundoff,  truncation,,  etc. 


2-  Analysis  of  Eesults 


Now,  the  influence  coefficients  are  actually  the 
various  partial  derivatives  of  the  argument  at  final  time 
with  respect  to  the  arguments  x,  y,  z,  u,  v,  and  w at 
initial  time.  For  instance,  the  influence  coefficient 
associated  with  x defines  a six  compbnent  vector  whose 
components  are  the  partial  of  x at  final  time  with  respect 
to  X at  initial  time,  the  partial  of  x at  final  time  with 
respect  to  y at  initial  time,  the  partial  of  x at  final  time 
with  respect  to  z at  initial  time,  the  partial  of  x at  final 
time  with  respect  to  u at  initial  tine,  the  partial  of  x at 
final  time  with  respect  to  v at  initial  time,  and  the 
partial  of  x at  final  time  with  respect  to  w at  initial 
tine.  The  values  of  all  these  partial  derivatives  or 
influence  coefficients  are  given  in  TABLE  V. 

These  'coefficients  serve  as  a cross  check  for  the 
analysis  of  the  variational  equations.  For  instance,  the 
variational  equations  with  DELu(0)=1.0  yielded  a value  of 
DELx  at  tf  of  32.867  while  the  adjoint  indicates  that  the 
partial  of  x at  final  time  with  respect  to  u at  initial  time 
is  32.870,  which  agree  to  four  figures.  Likewise,  the 
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variational  equations  with  OELv(0)-1.0  yielded  a DELy  at  t£ 
of  31.297  while  the  adjoint  indicates  that  the  partial  of  y 
at  final  time  with  respect  to  v at  initial  time  is  31.298. 
Other  values  were  checked  similarly. 

In  conclusion,  we  found  that  the  three  methods  of 
calculating  the  effects  of  Initial  conditions  on  the  point 
Qf  fall:  (1)  recomputing  the  trajectory,  (2)  using  the 
variational  equations,  and  (3)  using  the  adjoint  are 
satisfactory  within  the  range  of  interest.  iie  did  find  that 
the  differential  method  augmented  hy  the  correction  fcr  the 
change  in  terminal  time  of  fall  is  slightly  more  accurate, 
hut  the  increase  is  really  not  significant.  Finally,  we  see 
that  in  a computer  age  such  as  today  the  effects  can  just  as 
easily  be  calculated  by  changing  the  initial  condition  and 
programming  the  trajectory  again,  although  the  use  of  the 
variational  equations  is  still  an  important  tool  for 
analyzing  such  problems. 
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0.0 


0.0 


.920 

0.0 


0.0 

.667 


0.0 

0.0 

0.0 

o 

• 

o 

0.003 

.007 

o 

• 

o 

0.0 

.616 

,064 

. 147 

.301 

Table  V 

Adjoint  Sensitivity 

Factors 

33 

C.  COHMENTS 


The  pcoblems  of  Chapter  II  can  also  be  solved  by  the 
■ethods  of  this  chapter.  Indeed,  Bliss  developed  the 
techniques  associated  with  the  adjoint  to  handle  just  such 
problems,  while  he  was  at  Aberdeen  Proving  Grounds  during 
Uorld  War  I [Bef.  7]. 

Let  us  consider  the  effects  of  varying  the  wind  profile 
from  its  values  on  some  nominal  path.  The  variational 
equations,  (2a)  through  (2f) , apply  directly.  The  values  of 
z,  u,  Vwx,  etc.,  are  those  from  the  given  nominal  path,  and 
the  values  of  OELVwx,  OELVwy  are  the  values  on  the  new  path 
minus  those  cn  the  nominal  path. 

If  the  adjoint  equations  are  used,  the  equation  for,  say 
DELy(tf) , reduces  to 

where  the  solutions  to  the  adjoint  are  those  associated  with 
y,  for  which  P2(tf)=1,  and  the  final  value  of  the  other  Pi 
are  zero.  The  evaluation  of  the  integral  is  simplier  than 
integrating  a system  of  differential  equations. 
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IV.  ANALYSIS  QF  FACTORS  INFLOENC^G  AIRCRAFT  POSITION  MS 

VELOCITY  ESTIMATES 


The  effects  of  sose  errocsv  in  radar  aeasureoeots  of 
position  on  the  filtered  values  for  position  and  velocity 
are  investigated  here.  Also  included,  in  the  study  are 
errors  due  to  some  siaple  maneuvers  and  associated 
transients.  Reference  1 provides  the  equations  used  in  the 
AN/TPQ-27  system  to  predict  and  filter  the  aircraft  position 
and  velocity.  The  analysis  is  generally  carried  out  for  a 
single  coordinate,  y;  by  definition  xdot  is  zero  and 
ordinarily  zdot  will  also  be  zero. 


A.  RANDOM  UNCORRELATED  NOISE 


The  first  error  analyzed  is  that  due  to  random 
uncorrelated  errors  or  noise  in  the  radar  measurements.  The 
variance  of  the  resulting  errors  was  determined  in  terms  of 
the  variance  of  the  noise,  in  two  ways.  The  first  was 
purely  computational.  A set  of  random  numbers,  v (k) , with 
mean  value  zero,  veoriance  one,  and  gaussian  distribution  was 
generated.  These  were  used  like  measurements  or 
observations  in  the  filtering  equations.  Sequences  of 
fx.'tered  estimates  for  y and  ydot  were  generated  and  their 
means  (which  were  near  zero)  and  variances  were  calculated 
on  the  computer.  The  second  method  was  to  analyze  the 
filtering  equations  to  estimate  the  Variances  of  the 
filtered  estimates  of  y and  ydot. 
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1 . A Nuaerical  Approach 


A routine  for  generating  randon  nunbers  was  called 
froB  the  coaputer  library.  It  generated  random  numbers 
V ('<) , . . . .etc.  These  values  were  used  as  a forcing 
function,  y (k) , in  the  filtering  equations.  The  position 
and  velocity,  y and  ydot,  were  initially  set  to  zero,  and 
then  their  filtered  values  calculated  iteratively, 
numerical  estimates  were  then  made  for  the  standard 
deviations  of  the  filtered  values  of  y and  ydot  using  those 
calculated  values  as  the  sample  [Ref.  5].  It  was  then  found 
with  this  strictly  numerical  analysis  that  the  standard 
deviation  of  the  sample  of  filtered  values  of  y was 
approximately  .245  ^ and  the  standard  deviation  of  the 
sample  of  filtered  values  of  ydot  was  approximately  .118^  . 


2.  An  An^ytic  Approach 


The  prediction  equations  are: 

At 


A 

A 

A 

# 

y a. 

y 

y 

A 

A 

0 

V * 

0 

y 

The  filtering  equations  are: 

A A J 


o 


AtyA»-/ 


/ ‘ y 

where  y 


Cy-  y 

' (y  - L., 


J 


0) 

(zj 

(s) 

(*J 


is  the  estimate  of  y(tn)  made  at  time  tm;  the 
corresponding  estimate  for  ydot  is  y ; alpha=. 1, 

beta=,0052,  and  DT^,  125  [Ref.  1].  Now,  t^^ Substitution  of 
the  right  side  of  (1)  into  (3)  , and  the  substitution  of  the 
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right  sides  of  (1)  and  (2)  into  (4)  yields: 


} * f 


A / A O 

^ u At  ^ y ~ y 


-.(i-<^)t  +/it(i-<)  ^ 


%tLf-L  -y 


' -%t 


Y 


C/-  ^ y + 


/ + /ill*  y 

/H-/j/H~/  A>f  , 


These  are  the  equations  for  the  filtered  estimates  of  y and 
ydot,  in  terms  of  the  observations. 


Next,  the  following  substitutions  were  made: 


r 

^7j 

II 

. r 

)•  ^ 

{ ‘ y j 

'^j'n 

then  (5)  and  (6)  become 

/*»-/  y»I«/  A| 

('/«>; 

]r=  Tt  f ^ 

^ /n-/  **•"/  ^ 

00 
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One  of  the  variables,  say,  zeta  aay  be  elioinated  as 
zolloHs.  shift  the  indices  in  (10)  and  (11)  by  one,  to  get 


r) 

+ -f-  ^1/ 

Ot) 

]r  = 

/»»-/  L 

Oi) 

f w 

<•»-  z 

jf\~Z  '*»•/  • 

When  zeta  is  eiizinated  among  these  equations  the 

equation 

for  eta  becoaes 

7?  = - 

Of) 

^ "r  # 

/»!-  2 

• 

If  we  set 

A - 

os) 

6 --  - 0-^) 

Oi) 

c = ^ 

Oy) 

this  becoaes 

Of) 

^ /•»-/  ^-2 

/*7  # 

la  the  same 

vay  the  equation  for  zeta 

may  be 

obtained. 

-fit/  * 7)1/ 

Of) 

/If-/  /n-i 
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where 


£)»  -yS  . 0-°) 

In  order  to  estimate  the  variance  of  eta  we  need  to 
express  eta  in  terms  of  the  uncorreiated  observations.  To 
do  this  let  us  shift  indices  in  equation  (18) ; 

/!??  *8);  i-.c\/  * cy  (2-0 

and  eliminate  from  (18)  with  this; 

(P-) 

/•»  /*»-2  ^ 

•f  (/(of  f C)V  + /9C  \/ 

An-f  /h-2 

Bepeated  substitutions  of  this  type  lead  to  a series  which 
may  be  written 


?9  . c(\/h‘  c}y  c(  -h  /9C)I/ 

At 

+/>  cjy  y. 

* /*i-i 

or 


S 

•f 


1/  ^ 


1/ 

Z 


-f- 


(zi) 
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The  coefficients  A,  B,  and  E are  defined  by  tne 
recursion  formulas  and  initial  values  indicated  below. 


/(  . M +8 

a^J 

8 . 8/1 

/IW-/ 

£',/!<■(+/{  C 

<t7j 

^ /»*»-/ 

/f  * / 

Ctf) 

C 

B^O 

<z^) 

o 

£»c( 

Oo) 

o 


If  we  square  both  sides  and  consider  the  expected 
values  we  get 


(iO 


A similar  operation  for  zeta  leads  to  a recursion 

formula 


r.  /!  /?  o 

where  the  A's  and  B's  have  been  previously  defined  and 


He  then  ger 


<3i) 

03) 

0*) 
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3*  Coaparison  of  fiesalts  Obtained 


The  series  represented  in  equations  (31)  and  (34) 
were  suamed  iteratively.  Calculated  standard  deviations  for 
the  sample  of  filtered  values  of  y and  of  ydot  were  found  to 
be  .275  ^ and  .095  ^ respectively.  The  two  methods  of 
analysis  yielded  values  differing  by  approximately  eleven 
(11)  pe?  cent  for  the  standard  deviation  of  y and  by 
approximately  twenty  (20)  per  cent  for  the  standard 
deviation  of  ydot.  The  difference  of  twenty  (20)  per  cent 
is  somewhat  disappointing. 


B.  PEBIODIC  ESfiOBS 


iaeasurements  are  also  often  suoject  to  some  sort  of 
periodic  error,  or  noise.  This  was  approximated  by  a sin 
wave  and  the  effects  on  the  filtered  values  of  y and  ydot 
estimated.  Various  frequencies  were  considered  and  the 
resulting  error  in  position  and  velocity  was  determined  for 
each. 

1 . Analysis 

The  filtering  equations  (10)  and  (11)  were  of  the 

form: 

/p  » ci  y-  1/ 

S • Oi) 
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where 


c^v 


c » /“  « (se) 

The  errors  in  aeasurenents  were  expressed  in  complex 

form 

1/  ^ CsfJ 

where  theta  is  the  frequency.  ilultiple  values  of  pi, 
differing  by  factors  of  two,  were  used  for  theta,  and  the 
aaplitudes  of  the  resulting  steady  state  solution  were 
calculated. 


The  solutions  were  assumed  to  have  the  same  form. 


/*> 

^ .Ce 

^ 2 


(TO) 

On) 


and  these  were  substituted  into  equations  (35)  and  (36).  It 


is  convenient  to  shift 

the  index  by  one  in 

(35)  and  (36)  . 

a.?^  >* 

^ ^1/ 

(n) 

«7 

r • 

(n) 

/ni-i 

n? 

He  also  shift  the  index  in  (40)  and  (41)  by  one. 

fl  -C,e. 

tn*l  • 

Substituting  (44)  and  (45)  into  (42)  and  (43)  respectively 


on) 

(ff) 


yields 


CH) 
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Now,  (46)  can  be  written 

a = act’*’'  (f*) 

and  (47)  can  be  written 

^ i- ^■4e.^‘  . fffj 

Equations  (48)  and  (49)  are  two  equations  in  two 
unknowns.  Cl  and  C2.  The  solutions  are 

c,‘(f+g'^/(h^  sC) 

where 

/^=  «r  Cos  26  <,o(  CoS  6 t cos  6 iSi) 

Q n S//p  3L9  - e(  s/'o  6 t a, Si/f  6 C^J 

fi  s CL  <,  i-  C.OS  26^  ^CL'f'K, )coS  6 

S * ^ S//p  S//76 

and 

' O'*  f^yo  * ^0 

where 

^ CoS  26  - a.  Tf^  c,os  6 ^ o(  Cos  6 

s/^  - Ay<  s//;s  - si/)6  . 
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Then  the  amplitudes  of  the  sin  waves  are 


jc^l  ,]/P7fyp7P^  . 

2 . The  flesults 


The  magnitudes  of  Cl  and  C2  were  calculated  on  the 
computer  for  seven  different  periods,  and  the  results  are 
shown  in  TABLE  VI. 

These  results  show  that  when  the  errors  in  the 
signal  are  sinusoidal,  the  errors  in  the  filtered  values  are 
also  sinusoidal.  The  amplitudes  of  errors  in  the  filtered 
values  are  much  smaller  than  the  measured  values  when  the 
period  is  small.  For  large  periods  the  amplitudes  of  tne 
errors  in  the  filtered  values  of  y are  approximately  egual 
to  the  measured  values  of  y.  Periods  of  order  eight  (8)  to 
sixteen  (16)  seconds  appear  to  be  the  critical  periods, 
where  maximum  errors  in  the  filtered  values  of  y and  ydot 
occur. 

As  an  example  of  the  significance  of  a sinusoidal 
error  in  the  measured  value  of  ydot,  suppose  we  have  a 
siuLsoidal  error  with  an  amplitude  of  ten  (10)  feet  per 
second  and  a period  of  eight  (8)  seconds.  This  would  result 
in  an  error  in  the  filtered  value  of  velocity  of 
approximately  four  (e)  feet  per  second,  which  would  cause  a 
miss  distance  of  approximately  125  feet. 
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THETA 

PERIOD  (secs) 

l«i( 

^2  / A't 

IT 

.25 

.051 

.024 

Trn 

.5 

.073 

.032 

"r/A 

1 

.134 

.056 

▼/s 

2 

.265 

.112 

V/16 

4 

.537 

.224 

•*r/32 

3 

1.034 

.392 

ir/hA 

16 

1.239 

.36 

T/128 

32 

1 .091 

.192 

Table  VI 

Calculated  Amplitudes  for  Various  Periods 
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C.  EFFECTS  CF  TOBNS 


Two  types  of  oaneuvers  are  analyzed  here. 

First,  during  the  final  stage  the  plane  can  be  expected 
to  aaice  a sequence  of  corrections,  which  we  nay  expect  to  be 
snail,  sharp  turns.  Hence,  we  need  to  know  the  errors 
introduced  in  the  velocity  and  position  estinates  by  such 
naneuvers.  As  a typical  turn  we  considered  a plane  flying 
at  350  feet  per  second  turning  at  a rate  of  three  (3) 
degrees  per  second  for  one  second.  This  yields  values 


/f*  / /e5<p  (<u*) 

y C^)  ' /f-i/f  A J A-Sf  (ioi) 

yCAJ  - ?./s- 

tC*:)  • /e-A  J Cte^) 


nornal  to  the  initial  notion.  The  values  along  the  initial 
course  are  not  changed  significantly. 

The  values  of  y and  ydot  and  the  filtered  values  of  y 
and  ydot  were  calculated,  and  the  errors  checJced.  The 
naxiaun  errors  were  of  the  order  .0005  feet  for  y and  .0008 
feet  per  second  for  ydot,  which  are  negligible. 

Second,  we  considered  a plane  in  a constant  turning 
novenent.  This  is  an  unlihely  naneuver  for  an  aircraft  to 
be  perforning,  since  the  airplane  flies  on  a straight  line 
except  for  ninor  adjustnents  in  the  precision  guidance 
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stage.  For  this  reason  higher  errors  were  expected. 

He  calculated  the  errors  in  the  filtered  values  of  y and 
ydot,  and  the  corresponding  miss  distances  as  a result  of 
these  naneuvers.  It  was  found  that  errors  in  the  filtered 
values  of  y caused  a miss  distance  of  approximately  100 
feet,  while  errors  in  the  filtered  values  of  ydot  caused  a 
miss  distance  of  approximately  1100  feet.  This  is  a 
significant  miss  distance,  but  as  stated  earlier,  it  is  not 
surprising  since  a constant  turn  by  the  plane  should  only 
occur  in  the  preliminary  stages. 
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V.  COMCLU SIGNS  AND  COUMfiNTS 


Several  conclusions  and  some  comments  ace  presented  here 
as  a means  ot  summarizing  this  study. 

This  thesis  had  several  purposes,  focusing  on  the  desire 
of  the  author  to  study  mathematical  methods  and  techniques 
that  have  been  and/or  will  be  useful  when  applied  to  a 
problem  of  interest  to  the  Marine  Corps.  The  techniques  and 
methods  used  here  are  all  known,  but  there  is  no  place  where 
the  combination  of  numerical  methods,  ballistics, 
differential  techniques,  difference  (filtering)  equations, 
and  statistics  is  available;  they  generally  are  scattered 
throughout  various  fields. 


A.  C^NCLUSICMS 


Most  of  the  results  obtained  in  this  thesis  were,  in  a 
sense,  expected.  The  range  of  validity  of  linear 
approximations  to  changes  was  surprisingly  large,  however. 

In  our  very  limited  checking,  the  fraction  eight^tenths 
of  wind  at  altitude,  actually  given  to  us  as  part  of  the 
"lore"  was  found  to  be  near  the  best  choice  for  use  in  the 
absence  of  stratified  wind  data.  Its  validity,  however,  can 
be  checked  only  by  some  actual  testing  and  studying  of 
actual  wind  structures. 

He  also  found  that  errors  in  various  initial  conditions 
result  in  errors  in  the  impact  point  of  the  bomb  which  are 
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linear  over  a fairly  vide  range.  It  does  not  seen  to  aake 
■ucb  difference  which  technigue  we  use  to  estimate  these. 
Under  the  conditions  of  our  trajectories  an  error  or  change 
of  one  foot  per  second  in  horizontal  velocity  led  to  an 
error  or  change  in  pq'int  of  impact  of  about  32  feet,  roughly 
Independent  of  the  direction  of  the  error.  A value  of  zdot 
of  one  (1)  foot  per  second  at  release  led  to  an  increase  in 
range  of  about  2.26  feet.  The  miss  distances  were  found  to 
be  nearly  linear  over  a wide  range  of  values  for  the  initial 
errors. 

The  analysis  of  the  difference  equations  was  used  to 
show  how  a random  error  in  measurement  leads  directly  to  a 
random  error  in  the  point  of  impact.  From  the  standard 
deviation  of  the  radar  measurements  we  derived  the  standard 
deviation  of  the  estimates  for  y and  ydot  and  their 
correlation.  With  these  and  the  sensitivity  coefficients  we 
obtained  the  standard  deviation  of  the  associated  misses. 
Under  the  conditions  we  found  that  if  the  radar  has  a random 
error  in  measurement  with  standard  deviation  ^ , then  there 
will  be  a resulting  error  in  point  of  fall  with  standard 
deviation^  example,  if  the  radar  noise  has  an 
rms  (root-mean-square)  value  of  ten  (10)  feet,  the  point  of 
impact  will  have  an  associated  rms  value  of  thirty  (30) 
feet,  dost  of  this  results  from  the  error  in  predicting 
ydot,  though  y and  ydot  have  a fciirly  large  positive 
correlation,  and  hence  the  error  is  larger  than  if  they  were 
independent. 

Also  we  found  that  the  presence  of  a sinusoidal  error  in 
the  measurements  resulted  in  errors  which  are  sin  waves  of 
varying  amplitude  in  the  position  and  velocity  estimates. 
He  found  that  a sinusoidal  error  with  an  amplitude  of  ten 
(10)  feet  and  a period  of  eight  (8)  seconds  will  cause  a 
periodic  error  in  the  point  of  fall  with  a maximum  amplitude 
of  125  feet  and  an  rms  value  of  about  96  feet.  The  actual 
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error  depends  on  the  particular  time  of  release.  Sinusoidal 
errors  with  periods  of  one  (1)  second  or  less  have  little 
effect. 

The  effects  of  a small,  sharp  turn  by  the  aircraft  on 
the  position  and  velocity  estimates  were  found  to  be  fairly 
negligible.  On  the  other  hand,  the  errors  in  position  and 
velocity  estimates  for  an  aircraft  in  a constant  turn  were 
found  to  be  so  large  that  they  would  generally  preclude  the 
use  of  the  present  filtering  equations  (For  a turn  rate  of 
three  (3)  degrees  per  second  a miss  distance  of  1100  feet 
would  result,  principally  from  the  error  in  the  filtered 
value  of  ydot.):  however,  this  type  of  maneuver  should 
never  occur  during  the  final  precision  guidance  stage. 


a.  COailENIS 


A number  of  comments  pertaining  to  our  study  should  be 
made.  First,  a number  of  errors  and  their  sources  were  not 
studied,  primarily  since  little  could  be  done  at  this  level 
from  a mathematical  or  a scientific  nature  to  check  or 
improve  those  errors,  such  errors  include  radar  alignment 
error,  the  actual  bomb  alignment  on  the  aircraft,  the 
dampness  or  dryness  of  the  bomb  itself,  and  the  behavior  of 
air  currents  in  the  atmosphere. 

Proper  radar  alignment  depends  upon  proper  training  and 
discipline  of  the  personnel  involved  in  setting  up  the 
equipment,  as  well  as  good  surveys  prior  to  setting  up  the 
radar.  These  are  factors  that  must  be  considered  in  the 
field  where  the  equipment  is  operated. 

The  assumption  in  this  thesis  that  the  tail  vanes  kept 
the  bomb  aligned  with  its  axis  tangential  to  the  velocity 
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vector  relative  to  the  air  is  felt  to  be  probably  the 
greatest  source  of  random  error  in  the  systea.  The  bomb 
cones  out  with  angular  notion  and  nay  even  tumble 
nonentarily;  while  this  occurs  at  the  tine  of  low  speed  and 
thinnest  air  the  velocity  changes  at  this  tine  cause  large 
integrated  effects.  Any  analysis  of  this  attempting  to 
model  the  aerodynamic  forces  and  moments  due  to  various 
initial  configurations  would  be  of  doubtful  value.  The  only 
feasible  way  seems  to  be  numerical  experimentation  with  the 
actual  configuration,  combined  with  dropping  the  bombs  and 
observing  the  paths. 

He  also  limited  ourselves  to  -one  set  of  standard 
conditions;  there  is  little  point  to  carrying  out  more 
computations  of  the  same  type  unless  the  computations  are 
associated  with  an  actual  system.  Then  many  of  these  should 
be  carried  out  over  the  range  of  parameters  expected  to 
occur,  such  as  altitude  of  plane  and  of  target;  speed  of  the 
plane,  and  particularly  with  the  best  approximation  for  the 
drag  coefficient,  as  a function  of  :iach  number. 

It  must  be  noted  that  this  research  took  a general  look 
at  a sophisticated  problem,  and  in  doing  so  touched  only  the 
surface  in  terms  of  investigating  error  sources.  A system 
as  expensive  and  as  sophisticated  as  the  AN/IPQ-27  requires 
intensive  research  and  testing,  as  it  has  received  and  will 
continue  to  receive. 

Finally,  an  Important  facet  of  this  work  was  to  provide 
a familiarization,  or  an  awareness  document  for  ready 
reference  to  those  individuals  who  need  some  introduction  to 
the  mathematics  used  to  analyze  the  system.  In  that  sense 
the  intended  goal  of  the  research  was  achieved. 
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